Spin Polarization of Quasi Two-Dimensional Hole Systems 
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In quasi two-dimensional (2D) hole systems with an effective spin j = 3/2, heavy hole-light hole 
splitting results in a quantization of angular momentum perpendicular to the 2D plane. The spin 
polarization of quasi 2D hole systems due to an in-plane magnetic field Bit thus competes with the 
heavy hole-light hole splitting. As a result, the spin polarization of hole systems is very different 
from the spin polarization of j — 1/2 electron systems. In particular, it is shown that the spin 
polarization of quasi 2D heavy hole systems can change its sign at a finite value of B« . 
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A spin splitting of the energy levels of the electrons in 
a solid can be caused by an external magnetic field B, the 
Zeeman splitting, and we can also have a spin splitting 
due to the inversion asymmetry (IA)of the crystal struc- 
ture or the confining potential [jj. The IA- induced 
spin splitting can be ascribed to an effective magnetic 
field |3j • Intimately related with the spin splitting of the 
energy levels is the spin orientation or spin polarization 
of the wave functions. These quantities are of great im- 
portance for our understanding of the fundamental prop- 
erties of electron systems. Recently, the interest in this 
subject has been renewed because it has been noted that 
the spin polarization of quasi two-dimensional (2D) sys- 
tems due to a magnetic field B» parallel to the 2D plane 
affects the apparent metallic behavior of these systems 
0. Furthermore, the spin polarization is an important 
parameter for possible applications in the field of spin- 
tronics [jj. We distinguish here between the spin ori- 
entation which is a property of the individual (occupied 
or empty) states, and the spin polarization which is the 
average spin orientation of all occupied states. In quasi 
2D systems the spin polarization has been studied by 
applying a ma gnet ic field parallel to the 2D plane 
0, El S H 01H El If the Zeeman en- 
ergy splitting due to Bn becomes sufficiently large, the 
minority spin subband is completely depopulated. In 
such a situation, electron systems are fully spin-polarized 

111100. 

Recently, there has been considerable interest in the 
spin polarization of quasi 2D hole systems 0,0,0,0]- 
Hole systems in the uppermost valence band of many 
common semiconductors such as Ge and GaAs are differ- 
ent from electron systems due to the fact that the hole 
states have an effective spin j = 3/2 0. Subband quan- 
tization in 2D hole systems yields a quantization of angu- 
lar momentum with z component of angular momentum 
m = ±3/2 for the so-called heavy hole (HH) states and 
to = ±1/2 for the light hole (LH) states. Often it is 
assumed that the spin degree of freedom of HH and LH 
states in quasi 2D hole systems behaves similar to the 
spin of j = 1/2 electron states. This assumption ap- 
pears to be justified when the HH-LH splitting is large. 
However, we will show here that such an approach can 



be misleading. The important point is that the quan- 
tization axis of angular momentum that is enforced by 
HH-LH splitting is pointing perpendicular to the plane 
of the quasi 2D system. On the other hand, in general 
the effective Hamiltonians for IA-induced spin splitting 
and Zeeman splitting at Bu > tend to orient the spin 
vector parallel to the plane of the quasi 2D system |3|. 
The energy scale for IA-induced spin splitting and Zee- 
man splitting are typically much smaller than the HH-LH 
splitting. Thus it has been noted that for HH states the 
IA-induced spin splitting E3 an d the Zeeman splitting 
at Bh > [l8| are suppressed, i.e., in a model with axial 
symmetry, the IA-induced spin splitting and the Zeeman 
splitting are higher-order effects proportional to fc| and 
Bu, respectively EH- Here fey is the in-plane wave vec- 
tor. This result reflects the fact that we can have only 
one quantization axis of angular momentum. A detailed 
discussion of HH-LH splitting can be found in Ref. [2fJ • 

In the present work we thus want to investigate 
the spin orientation and spin polarization that can be 
achieved in quasi 2D hole systems. We will show that 
in quasi 2D hole systems these quantities behave funda- 
mentally different from what is well-known about these 
quantities in spin-1/2 electron systems. Throughout, we 
will assume that only one orbital subband is occupied. 

First, we briefly review for comparison the main fea- 
tures of spin orientation and spin polarization in electron 
systems. In j = 1/2 electron systems at B = we have 
two degenerate states |f) and ||) that correspond to spin 
up and spin down. A finite magnetic field B gives rise to a 
Zeeman term cx er • B in the Hamiltonian which results in 
a spin splitting cx ±B of the energy levels. Here B can be 
either an external magnetic field or the k-dependent ef- 
fective magnetic field due to a spin-orbit interaction such 
as the Dresselhaus El or Rashba term 0. The symbol 
<r denotes the vector of Pauli spin matrices. The eigen- 
states |±) of the Hamiltonian at B > are orthonormal 
linear combinations of the states {||) , ||)} with spin ori- 
entation S + = (+|er|+) parallel to B and S_ = (— \tr\— ) 
antiparallel to B. Furthermore, it is easy to show that for 
normalized states |±) the spin orientations S± are always 
unit vectors Q. We see here that the spin orientation 
vectors S± of electron states always follow the magnetic 
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field B. The spin polarization at B > is obtained by 
averaging over the spin orientation of all occupied states. 
In the presence of an external magnetic field By parallel 
to the 2D plane, the system will be fully spin polarized 
if the Zeeman energy splitting is large enough to com- 
pletely depopulate the minority spin subband If By 
is an effective magnetic field describing the spin splitting 
in inversion-asymmetric systems, the k|| dependence of 
By is such that the spin polarization is zero ||. 

Hole states with an effective spin j — 3/2 are charac- 
terized by the 4x4 Luttinger Hamiltonian 0] . The cor- 
responding wave functions are four-component spinors. 
In quasi 2D systems we get HH and LH states as dis- 
cussed above. At B = the HH and LH states are each 
two-fold degenerate. We can choose the eigenstates 



where 



|* L 5 



) = 



2tt 



2tt 



/ 1 \ 



«4 



1 






1 

±1 
V o 1 



( \ 



(la) 



(lb) 



Here we have assumed that the Luttinger Hamiltonian 
is expressed in a basis of j = 3/2 angular momentum 
eigenf unctions in the order m = 3/2, 1/2, —1/2, and 
—3/2; and the quantization axis of these basis functions 
is perpendicular to the 2D plane. In Eq. the sym- 
bol a denotes the subband index. For simplicity we have 
neglected HH-LH mixing due to a nonzero in-plane wave 
vector k|| and we also neglected the k|| dependence of the 
envelope functions £ a (z). These aspects are fully taken 
into account in the numerical calculations discussed be- 
low. One can easily convince oneself that, in contrast 
to j = 1/2 electron systems, one can only find a lin- 
ear combination of the degenerate HH states 1^™) and 
that are eigenstates of J z , but there are no such 
linear combinations that represent eigenstates of J x or 
J y . Here J x , J y , and J z are angular momentum matri- 
ces for j — 3/2. This fact clearly illustrates the frozen 
angular momentum of HH states. 

A finite in-plane magnetic field Bu gives rise to a mix- 
ing of HH and LH states. In first order degenerate per- 
turbation theory we obtain the following expressions for 
the eigenstates of quasi 2D holes in an in-plane field Bu 
in x direction 
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are positive constants, k is the isotropic g factor and 
72 is the second Luttinger parameter in the Luttinger 
Hamiltonian /is is the Bohr magneton, m is the 
mass of free electrons, E% (E l a ) is the energy of the cuth 
HH (LH) subband at Bii — 0, and we have used the 
gauge A = (0, — zB\\ , 0). In our notation, the HH "+" 
states (a = 1) belong to the majority spin subband that 
remains occupied for large By . For simplicity, we have 
neglected the small terms proportional to the anisotropic 
g factor q. We also neglected the coupling to neighboring 
subbands a' ^ a. We note that JC oc w 2 and Q cx w 4 , 
where w is the width of the quasi-2D layer. 

Similar to electron states the spin orientation of hole 
states is given by 



S Q± = 3(* Q± |J|* Q± ) 



(4) 



Here J = (J x , Jy, Jz) is the vector of angular momen- 
tum matrices for j = 3/2 and the prefactor 2/3 ensures 
that | S a ±| < 1- F° r a magnetic field in x direction only 
the x components of the vectors S Q ± are nonzero. (As 
discussed below we also get a y component if the k|| de- 
pendence of \^ a ±) is taken into account.) Using Eq. (J2J 
we obtain in first order of the HH-LH splitting E% — E a 
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For typical values of Bii the second and third term in Eq. 
(|5b|) are much smaller than 2/3. Thus, similar to elec- 
tron states at Bu > 0, the spin orientations and 5^5 
of the Zeeman-split LH states are essentially antiparallel 
to each other and independent of B\\. For small fields 
B\\ the Zeeman-split HH states have the same spin ori- 
entation antiparallel to By which is "compensated" by a 
term of opposite sign in Eq. (|5 b|) . The spin orientations 
S 1 ™ vanish for Bu — > 0. We see here that HH states 
are fundamentally different from j = 1/2 electron states 
where the spin orientation vectors of the spin-split states 
are always antiparallel to each other, and they have al- 
ways length one, independent of the magnitude of Bu 0. 
According to Eq. (|5a|l the spin orientation vectors of HH 
states are the smaller in magnitude the larger the HH-LH 
splitting E% — E l a , which is opposite to what one expects 
according to the picture that the spin of hole systems 
with a large HH-LH splitting behaved similar to the spin 
of j = 1/2 systems. Once again, these results reflect the 
frozen angular momentum of hole states due to the HH- 
LH splitting. For larger fields Bu the term Q quadratic 
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in Bn dominates over the linear term JC in Eq. J5a}. For 
the HH "+" states, we thus obtain a sign reversal of the 
spin orientation 5™ at a finite value of B\\ . The term K, 
is proportional to the hole g factor k whereas Q is pro- 
portional to the Luttinger parameter 72 characterizing 
the orbital motion of the holes. We can thus attribute 
the sign reversal of S 1 ™ to a competition between the 
effect of the motion of the orbital and the spin degree 
of freedom. Wc remark here that usually only the first 
HH subband is occupied in 2D hole systems. The sign 
reversal of S 1 ™ as a function of Bn is thus an important 
feature of these systems. It depends sensitively on the 
width w of the quasi-2D layer. 

The functional form of the Zeeman energy splitting 
in an in-plane magnetic field Bu is very different from 
the expressions for the spin orientation in Eq. ©. We 
emphasize here that the spin orientation depends on the 
eigenf unctions of the Schrodinger equation whereas the 
Zeeman splitting depends on the eigenvalues. In second 
order perturbation theory the Zeeman energy splitting 
due to an in-plane field B11 in x direction reads [2fJ 
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The product of n and 72 in Eq. JJJ as well as the cubic 
dependence on Bu indicate that in contrast to the spin 
orientation of HH states the Zeeman energy splitting is 
due to the combined effect of the motion of the orbital 
and the spin degree of freedom. The expressions for the 
energy eigenvalues also contain terms of even order in Bu . 
However, these terms do not contribute to the Zeeman 
splitting. We note that Z oc w 4 . 

If the spin splitting is caused by the ky-dependent ef- 
fective magnetic field in inversion- asymmetric systems we 
obtain results similar to Eqs. (J5J and © with Bn replaced 
by fcii . In particular, once again the spin vectors S a ± (kii ) 
of the spin-split HH states for a given in-plane wave vec- 
tor k|| are in general neither unit vectors nor they are 
antiparallel to each other. 

We remark that previously it has proven useful to dis- 
cuss spin splitting in quasi 2D HH systems by means of 
2x2 Hamiltonians acting only in the subspace of HH 
states 20]. On the other hand, in the present discussion 
it is crucial that we use the full 4x4 Luttinger Hamilto- 
nian to evaluate the spin orientation of quasi 2D systems. 
Here it should be kept in mind that the 2x2 Hamiltonians 
are effective Hamiltonians that emerge from a perturba- 
tive unitary transformation of the full 4x4 Luttinger 
Hamiltonian in order to decouple the motion of HH and 
LH states. To evaluate the spin orientation S we must 
either express J in the transformed basis or we must 
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FIG. 1: Spin orientation as a function of kii for a sym- 
metric (100) GaAs-Alo.3Gao.7As quantum well with hole den- 
sity p = 1 x 10 11 cm -2 and different values of an in-plane mag- 
netic field Bu in x direction. The left (right) column refers 
to the majority (minority) spin subband. The dimensions of 
the arrows are proportional to S™(ky). The dashed lines 
indicate the Bu -dependent Fermi contours. 



(re-)express the HH wave functions as four-component 
spinors as in Eqs. and (J5J. 

To confirm the qualitative results in Eq. (J5J we show 
in Fig. ^ the self-consistently calculated spin orienta- 
tion S 1± as a function of k|| for a symmetric (100) 
GaAs-Alo.3Gao.7As quantum well with hole density p = 



1x10 



M 2 and different values of an in-plane magnetic 
The left (right) column refers 
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FIG. 2: (a) Spin subband density p_ of the HH minor- 
ity spin subband and (b) spin polarization 7r as a func- 
tion of the in-plane field B\\ calculated for a symmetric 
(100) GaAs-Alo.3Gao.7As quantum well with hole density 
p = lx 10 11 cm -2 . 

to the majority (minority) spin subband. The dashed 
lines show the Fermi contours. The numerical calcula- 
tions follow Ref. (23. They fully take into account the 
anisotropy of the crystal structure. We see in Fig.^that 
the magnitude and the orientation of S 1 J t (kii) depend 
sensitively on Bu and k|| , in sharp contrast to electron 
states. For small values of Bu and fcii , the spin orien- 
tation vectors of the HH states in both spin subbands 
are essentially antiparallel to B||, in agreement with Eq. 
(|5a|l . For larger values of fc|| we can have a finite angle 
between S™(k||) and B||. Furthermore, the spin vectors 
S™(k||) and S™(k||) for a given wave vector k|| are in 
general neither parallel nor antiparallel. We remark that 
the k|| dependence of the spin orientation can be studied 
analytically using the model underlying Eqs. (0) and JSJ). 
However, we do not reproduce here the lengthy formulas. 

For the same system as in Fig. ^ we show in Fig. [2Ja) 
the self-consistently calculated depopulation of the HH 
minority spin subband, and Fig. 13b) shows the x com- 
ponent of the normalized spin polarization 

7T = - f dk l e \ E -± ( k ll ) - M S q± (k,| ) . (8) 

P a,± 

Here E a ±(k.n) is the spin-split subband dispersion and 
E-p is the Fermi energy. If B|| is pointing in x direc- 
tion, the y and z components of 7r are zero. In Eq. JSJ) 
we have assumed that hole energies are negative and the 
temperature is zero. We evaluate Eq. © by means of 
analytic quadratic Brillouin zone integration |2l|. It be- 
comes clear from Fig. |2 that the depopulation of the HH 



minority spin subband due to _B|| does not imply full spin 
polarization of the system. In Fig. [5] the minority spin 
subband is completely depopulated at Bo ss 28.2 T while 
tt(Bd) ~ 0.14. In agreement with Eq. I|5a|) we obtain a 
sign reversal of it at B\\ w 21.4 T. The derivative of 7r(B|| ) 
is discontinuous at Bu = Br). 

We have checked that the surprising results reported 
here are observed for a large variety of 2D HH systems, al- 
though the details depend on system parameters such as 
the density p, the width and the crystallographic growth 
direction of the quasi 2D system, and the material com- 
position of the layers. In the calculations, we have ne- 
glected many-particle effects which will be important for 
a quantitative analysis of spin polarization in the density 
regime discussed here. However, to the best of our knowl- 
edge an appropriate theory for particles with j — 3/2 is 
presently not available. We expect that the essential fea- 
tures in Figs. ^ and |2 will be confirmed by such a more 
quantitative theory. 
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